LBNL-41826 
hep-ph/9805405 

Radiative Decay of a Long-Lived Particle and Big-Bang 

Nucleosynthesis 



Erich Holtmann'*, Masahiro Kawasaki^, Kazunori Kohri*, and Takeo Moroi'^ 
Theoretical Physics Group, Lawrence Berkeley National Laboratory, 
University of California, Berkeley, California 94720 
Institute for Cosmic Ray Research, The University of Tokyo, Tanashi 188-8502, Japan 



Abstract 

The effects of radiatively decaying, long-lived particles on big-bang nucleosyn- 
thesis (BBN) are discussed. If high-energy photons are emitted after BBN, 
they may change the abundances of the light elements through photodisso- 
ciation processes, which may result in a significant discrepancy between the 
BBN theory and observation. We calculate the abundances of the light el- 
ements, including the effects of photodissociation induced by a radiatively 
decaying particle, but neglecting the hadronic branching ratio. Using these 
calculated abundances, we derive a constraint on such particles by compar- 
ing our theoretical results with observations. Taking into account the recent 
controversies regarding the observations of the light-element abundances, we 
derive constraints for various combinations of the measurements. We also 
discuss several models which predict such radiatively decaying particles, and 
we derive constraints on such models. 
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I. INTRODUCTION 



Big-bang nucleosynthesis (BBN) has been used to impose constraints on neutrinos and 
other hypothetical particles predicted by particle physics, because BBN is very sensitive to 
the thermal history of the early universe at temperatures T < 1 MeV [1]. 

Weakly interacting, massive particles appear often in particle physics. In this paper, 
we consider particles which have masses of ~ 0(100 GcV) and which interact with other 
particle only very weakly {e.g., through gravitation). These particles have lifetimes so long 
that they decay after the BBN of the light elements (D, ^He, '^He, etc.), so they and their 
decay products may affect the thermal history of the universe. In particular, if the long-lived 
particles decay into photons, then the emitted high-energy photons induce electromagnetic 
cascades and produce many soft photons. If the energy of these photons exceeds the binding 
energies of the light nuclides, then photodissociation may profoundly alter the light element 
abundances. Thus, we can impose constraints on the abundance and hfetime of a long-hved 
particle species, by considering the photodissociation processes induced by its decay. There 
are many works on this subject, STich as the constraints on massive neutrinos and gravitinos 
obtained by the comparison between the theoretical predictions and observations [2-6].^ 

A couple of years ago, Hata et al. [8] claimed that light-element observations seemed to 
conflict with the theoretical predictions of standard BBN. Their point was that standard 
BBN predicts too much ^He, if the baryon number density is determined by the D abundance 
inferred from observations; cquivalently, standard BBN predicts too much D, if the baryon 
number density is determined by the '^He observations. Inspired by this "crisis in BBN," 
many people re-examined standard and non-standard BBN by including systematic errors 
in the observations, or by introducing some non-standard properties of neutrinos [9,10]. In 
a previous paper [11], we investigated the effect upon BBN of radiatively decaying, massive 
particles. These particles induce an electromagnetic cascade. We found that in a certain pa- 
rameter region, the photons in this cascade destroy only D, so that the predicted abundances 
of D, '^Hc, and "^He fit the observations. 

However, since the "BBN crisis" was claimed, the situation concerning the observations of 
deuterium has changed. The D abundances in highly red-shifted quasar absorption systems 
(QAS) have been observed. The abundance of D in high- 2; QAS is considered to be the 
primordial value. Thanks to these direct new observations, we no longer need to use poorly 
understood models of chemical evolution to infer the primordial abundance from the material 
in solar neighborhood. 

Moreover, there are also differing determinations of the primordial ^He abundance. Hata 
et al. used a relatively low ^He abundance (viz., Y ~ 0.234, where Y is the primordial 
mass fraction of *^He) [12,13]. However, a higher *^He abundance {Y ^ 0.244) has also 
been reported [14-16], and it has been noted that this higher observation alleviates the 



^As pointed out in Ref. [7], even if the parent particle decays only into photons, these photons 
will produce hadrons with a branching ratio of at least 1%. However, since there is no data 

on some crucial cross sections involving '''Li and ^Bc, wc cannot include hadrodissociation in our 
statistical analysis. Since we have neglected hadrodissociation, our constraints may be regarded as 
conservative bounds. 
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discrepancy with standard BBN theory [17]. The typical errors in ^He observations are less 
than ~ 0.005, so we have discordant data for ^He. 

Since we have discordant "^He abundances and new observations for D, the previous 
constraint on the radiative decay of long-lived particles must be revised. In addition, the 
statistical analyses on radiatively decaying particles are insufficient in the previous works. 
Therefore, in our present paper, we perform a better statistical analysis of long-lived, ra- 
diatively decaying particles, and of the resultant photodissociations, in order to constrain 
the abundances and lifetimes of long-lived particles. In deriving the constraint, we use both 
high and low values of the ^He abundance, because it is premature to decide which data 
are correct. As a result, it will be shown that for low values of the ^He abundance, we have 
a poor agreement between the observations and the standard BBN theory. Moreover, we 
show in this case that a long-lived particle with appropriate abundance and lifetime can 
solve the discrepancy. In the case of high "^He, standard BBN fits the observations, so we 
derive stringent constraints on the properties of long-lived particles. 

In this paper, we also include the photodissociations of ^Li and ^Li for the first time. 
As we will show later, the destruction of ^Li docs not dramatically affect the predicted D 
and ^He, in the region where the observed D and '^He values are best fit. However, the ^Li 
produced by the destruction of ^Li can be two orders of magnitude more abundant than the 
standard BBN prediction of ^Li/H ~ O(10~^^). We discuss the possibihty that this process 
may be the origin of the ^Li which is observed in some low-metallicity halo stars. 

In Sec. II, we study how consistent the theoretically predicted abundances and obser- 
vations are, in the case of standard BBN. The radiative decay of long-lived particles is 
considered in Sec. Ill, and the particle physics models which predict such long-lived parti- 
cles are presented in Sec. IV. Finally, Sec. V is devoted to discussion and the conclusion. 



II. STANDARD BIG-BANG NUCLEOSYNTHESIS 

We begin by reviewing standard big-bang nucleosynthesis (SBBN). We are interested in 
the light elements, since their primordial abundances can be estimated from observations. In 
particular, we check the consistency between the theoretical predictions and the observations 
for the following quantities: 

^2 = rio/riH, (1) 

Y = p4He/pB, (2) 

ye^neu/nn, (3) 
r = mHe/riD, (4) 
Vv = nru/nji, (5) 

where pe is the total baryon energy density. 

In this section, we first review the observations of the fight elements, and the extrapola- 
tions back to the primordial abundances. Next, we describe our theoretical calculations of 
these abundances, by using standard big-bang theory as an example. Finally, we compare 
the theoretical and observed light-element abundances to determine how well the SBBN 
theory works. 
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A. Review of Observation 



Let us start with a review of the observations of the hght-element abundances. Two 
factors comphcate the interpretation of the observations of the hght-element abundances. 
First, there are various observational determinations for ^He which are not consistent with 
each other, within the quoted errors. This fact suggests that some groups have underesti- 
mated their systematic error. ^ We believe it is premature to judge which measurements are 
reliable; hence, we consider both of the observations when we test the consistency between 
theory and observation. Second, some guesswork is involved in the extrapolation back from 
the observed values to the primordial values, as we shall discuss below. Keeping these factors 
in mind, we review the estimations of the primordial abundances of D, ^He, ^He, ^Li, and 
^Li. 

D/H has been measured in the absorption lines of highly red-shifted (and therefore 
presumably primordial) Hi (neutral hydrogen) clouds which are backlit by quasars. The 
latest result suggests [19] 

" = (3.39 ± 0.25) X 10"^ (6) 

We use this value rather than the higher abundance which had been reported [20-23] in 
some of the old measurements of D/H. The older results suggested the abundance was 
y2 ~ O(10~^). However, we believe these results to be much more uncertain. For example, 
the authors of Ref. [22] admitted a large uncertainty in their results. Furthermore, results 
given in Ref. [23] are based on the fit of only the Lyman alpha limit, and the resolution is 
not good. Therefore, we will not use the high D values in deriving the constraints, but we 
will just discuss the implications of taking the high value seriously. 

For ^He, we use the pre-solar measurements. In this paper, we do not rely upon any 
models of galactic and stellar chemical evolution, because of the large uncertainty involved 
in extrapolating back to the primordial abundance. But it is reasonable to assume that 
^He/D is an increasing function of time, because D is the most fragile isotope, and it is 
certainly destroyed whenever ^He is destroyed. Using the solar-system data reanalyzed by 
Geiss [25], the ^He/D ratio is estimated to be [26] 

t''^^ ^ ivt/yt) ^ = 0.591 ± 0.536, (7) 

where denotes the pre-solar abundance. We take this to be an upper bound on the 
primordial ^He to D ratio: 

r"^' < r°^'. (8) 

Because the theoretical prediction of ^He/D in SBBN agrees so well with this upper bound, 
we do not include this constraint in the SBBN analysis. But when we investigate the 
photodissociation scenario, the situation is quite different. ^He photodissociation produces 



^It is also possible that primordial nucleosynthesis was truly inhomogeneous [18]. However, in 
this paper we adopt the conventional belief that BBN was homogeneous. 
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both D and ^He and can raise the ^He to D ratio [26]. Hence, in our analysis of BBN with 
photodissociation, we inchide this upper bound, as described in the appendix (sec Eq. (A5)). 
An analysis based upon the chemical evolution of ^He and D will appear in a separate paper 
by one of the authors [24] . 

The primordial ^He abundance is deduced from observations of extragalactic Hn regions 
(clouds of ionized hydrogen). Currently, there are two classes of Y°^^, reported by several 
independent groups of observers. Hence, we consider two cases: one low, and one high. 

We take our low "^He abundance from Olive, Skillman, and Steigman [13]. They used 
measurements of ^He and 0/H in 62 extragalactic Hn regions, and linearly extrapolated 
back to 0/H= to deduce the primordial value 

Low: Y"'' = 0.234 ± (0.002)3^* ± (O.OOS),^,.*. (9) 

(When they restrict their data set to only the lowest metallicity data, they obtain Y°'^^ = 
0.230 ± 0.003.) Their systematic error comes from numerous sources, but they claim that 
no source is expected to be much more than 2%. In particular, they estimate that stellar 
absorption is of order 1% or less. 

We take our high ^He abundance from Thuan and Izotov [15]. They used measurements 
of ^He and 0/H in a new sample of 45 blue compact dwarf galaxies to obtain 

High: = 0.244 ± (0.002),^^ ± (0.005),^,,^. (10) 

The last error is an estimate of the systematic error, taken from Izotov, Thuan, and Lipovet- 
sky [16]. Thuan and Izotov claim that Hej stellar absorption is an important effect; this ex- 
plains some of the difference between their result and that of Olive, Skillman, and Steigman. 

Rather than attempting to judge which group has done a better job of choosing their 
sample and correcting for systematic errors, we prefer to remain open-minded. Hence, we 
shall use both the high and low ^He abundances, without expressing a preference for one 
over the other. 

The ''Li/H abundance is taken from observations of the surfaces of Pop II (old generation) 
halo stars. ^Li is a fragile isotope and is easily destroyed in the warmer interior layers of a 
star. Since more massive (or equivalently, hotter) stars are mixed less, one might hope that 
the surfaces of old, hot stars consist of primordial material. Indeed, Spite and Spite [27] 
discovered a "plateau" in the graph of ''Li abundance vs. temperature of old halo stars, 
at high temperature. This plateau is interpreted as evidence that truly primordial ^Li has 
been detected. Using data from 41 plateau stars, Bonifacio and Molaro [28] determine the 
primordial value logio(?/7''^) — —9.762 ± (0.012) stat ± (0.05) syst- Bonifacio and Molaro argue 
that the data provides no evidence for ^Li/H depletion in the stellar atmospheres (caused 
by, e.g., stellar winds, rotational mixing, or diffusion). However, for our analysis, we shall 
adopt the more cautious estimate of Hogan [29] that ^Li may have been supplemented (by 
production in cosmic-ray interactions) or depleted (in stars) by a factor of two: [30] 

loglo(yf = -9-76 ± {0.012)stat ± {0.05)syst ± (O-S) factor of 2- (H) 

Because ^Li is so much rarer than ^Li, it is much more difficult to observe. Currently, 
there is insufficient data to find the "Spite plateau" of ^Li. However, we can set an upper 
bound on ^Li/^Li, since it is generally agreed that the evolution of ^Li is dominated by 
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production by spallation (reactions of cosmic rays with the interstellar medium). The upper 
bounds on ^Li/^Li observed in low-metallicity ([Fe/H] < —2.0) halo stars range from [31] 
y&lvi ~ 0.045 to ye/^/T ~ 0.13. (Note that the primordial ^Li and ^Li have both been 
destroyed in material which has been processed by stars.) 

Rotational mixing models [32] yield a survival factor for ^Li of order 0.05 and a survival 
factor for ^Li of order 0.005. Therefore, the upper bound for primordial ^Li/^Li ranges 
approximately from 

yt/yf^<Q.b to 1.3. (12) 

Since we have only a rough range of upper bounds on ^Li, and no lower bound, we will 
not use ^Li in our statistical analysis to test the concordance between observation and 
theory. Instead, we will just check the consistency of our theoretical results with the above 
constraint. 



B. Theoretical Calculations 

Theoretically, the primordial abundances of the light elements in SBBN depend only upon 
a single parameter: the baryon-to-photon ratio r}. In our analysis, we modified Kawano's 
nucleosynthesis code [33] to calculate the primordial light-element abundances and uncer- 
tainties. 

In our calculation, we included the uncertainty in the neutron lifetime [34], in the rates 
of the 11 most important nucleosynthesis reactions [35], and in the rates of the photofis- 
sion reactions (see Table II). We treated the neutron lifetime, the nucleosynthesis reaction 
rates, and the photofission reaction rates as independent random variables with Gaussian 
probability density functions (p.d.f.'s). We performed a Monte- Car lo^ over the neutron hfe- 
time, nucleosynthesis reaction rates, and photofission reaction rates, and we found that the 
light-element abundances were distributed approximately according to independent, Gaus- 
sian p.d.f.'s. Therefore, the p.d.f. p^^^ for the theoretical abundances is given by the product 
of the Gaussian p.d.f.'s 
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(13) 



for the individual abundances: 



th 



,th\ 



pGauss(^yth^ ^pGauss^Y'^^ X p^""^'* (log^Q 2/f ) • 



(14) 



In Fig. 1, we have plotted the theoretical predictions for the light-element abundances 
(solid lines) with their one-sigma errors (dashed lines), as functions of r]. 

The dependences of the abundances on rj can be seen intuitively [1,37]. The "^He abun- 
dance is a gentle, monotonically increasing function of rj. As rj increases, ^He is produced 



^It has recently been demonstrated that the uncertainties in SBBN can be quantified by the much 
quicker method of Unear propagation of errors. [36] 
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earlier because the "deuterium bottleneck" is overcome at a higher temperature due to the 
higher baryon density. Fewer neutrons have had time to decay, so more ^He is synthesized. 
Since ^He is the most tightly bound of the light nuclei, D and ^He are fused into "^He. The 
surviving abundances of D and ^He are determined by the competition between their de- 
struction rates and the expansion rate. The destruction rates are proportional to rj, so the 
larger r) is, the longer the destruction reactions continue. Therefore, D and ^He arc mono- 
tonically decreasing functions of rj. Moreover, the slope of D is steeper, because the binding 
energy of D is smaller than '^He. 

The graph of ^Li has a "trough" near rj ~ 3 x 10~^°. For a low baryon density 77 < 
3 X lO""^*^, ^Li is produced by ^He(T, 7)^Li and is destroyed by ^Li(p, a^Re. As rj increases, 
the destruction reaction become more efficient and the produced ^Li tends to decrease. On 
the other hand for a high baryon density ^ 3 x 10^^°, '^Li is mainly produced through 
the electron capture of ^Be, which is produced by ^B.e{a, 7)''Be. Because ''Be production 
becomes more effective as rj increases, the synthesized ^Li increases. The "trough" results 
from the overlap of these two components. The dominant source of ^Li in SBBN is D(q;, 
7)^Li. Thus, the rj dependence of ^Li resembles that of D. 

We have also plotted the 1-a observational constraints. The amount of overlap of the 
boxes is a rough measure of the consistency between theory and observations. We can also 
see the favored ranges of rj. However, we will discuss the details of our analysis more carefully 
in the following section. 



C. Statistical Analysis and Results 

Next, let us briefly explain how we quantify the consistency between theory and obser- 
vation. For this purpose, we define the variable 

(nth _ obs\2 

y2^V^-^ ' ' (15) 

where = (1/2, ^5 logio yr), and we add the systematic errors in quadrature: (0"°'"')^ = 
^^systy _|_ (^^staty^ j-ggg ^YiQ appendix for a detailed explanation of our use of x^.) ^ 
depends upon the parameters of our theory (viz. 77 in SBBN) through a*'* and cr*'*. 

Notice that we do not include ^Li in the calculation of x^, since the ^Li abundance has 
not been measured well. Instead, we check that y*^ /y\^ satisfies the bound (12). In the case 
of SBBN, we found that the ^Li abundance is small enough over the entire range of r] from 
8.0 X 10^^^ to 1.0 X 10~^. (Numerically, yQ^/yf' < 5 x 10^"^, which is well below the bound 
(12).) 

With this variable, we discuss how well the theoretical prediction agrees with observa- 
tion. More precisely, we calculate from the confidence level (CL.) with which the SBBN 
theory is excluded, at a given point in the parameter space of our theory (for three degrees 
of freedom) : 

C-L- = r l^HJ^y''^~''dy (16) 

(17) 
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In Fig. 2, we have plotted the and confidence level at which SBBN theory is excluded 
by the observations, as a function of rj. We find that high ^He is allowed at better than the 
68% C.L. at ?7 ~ 5 X 10~^°, while for low ''He, no value of rj works at the 91.5% C.L. 

The case of low ^He suggests a discrepancy with standard BBN. Some people beheve 
that this casts doubt on the low D or low ^He measurements [38]. However, we do not want 
to assume SBBN theory and use it to judge the validity of the observations; rather, we use 
the observations to test BBN theory. Therefore, we give equal consideration to both the 
high and low observed abundances of ^He. 

Before closing this section, we apply our analysis to constrain the number of neutrino 
species. Here, we vary r] and the number N^, of neutrino species, and we calculate the 
confidence level as a function of these variables. The results are shown in Fig. 3a,b. We can 
see that the standard scenario {Ny = 3) results in a good fit with a high value of ^He, while 
the case of low ''He results in a discrepancy. In fact, low ''He prefers N^, ~ 2, as pointed out 
by [8,10]. In Table I, we show the 95% C.L. bounds for the number of neutrino species N^, 
and T) in the two cases of the ^He abundances. 

III. BBN + X 

In this section, we discuss the implications of a radiatively decaying particle X for BBN. 
For this purpose, we first discuss the behavior of the photon spectrum induced by X. Then 
we show the abundances of the light elements, including the effects of the photodissociation 
induced by X. Comparing these abundances with observations, we constrain the parameter 
space for rj and X. 

A. Photon Spectrum 

In order to discuss the effect of high-energy photons on BBN, we need to know the shape 
of the photon spectrum induced by the primary high-energy photons from X decay. 

In the background thermal bath (which, in our case, is a mixture of photons 7bg, electrons 
CgQ, and nucleons Nbg), high-energy photons lose their energy by various cascade processes. 
In the cascade, the photon spectrum is induced, as discussed in various literature [39] . The 
important processes in our case are: 

• Double-photon pair creation (7 -|- 7bg — > e"^ -|- e~) 

• Photon-photon scattering (7 -|- 7bg — 7 + 7) 

• Pair creation in nuclei (7 -|- A'^bg + e~ + N) 

• Compton scattering (7 -|- e^Q — > 7 -|- e~) 

• Inverse Compton scattering (e^ -|- 7bg ^ -|- 7) 

(We may neglect double Compton scattering 7 + e^Q ^ 7 + 7 + 6", because Compton 
scattering is more important for thermalizing high-energy photons.) In our analysis, we 
numerically solved the Boltzmann equation including the above effects, and obtained the 
distribution function of photons, f^{E^). (For details, see Refs. [4,5].) 
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In Fig. 4, we show the photon spectrum for several temperatures T. Roughly speaking, 
we can see a large drop-off at ~ mfj22T for each temperature. Above this threshold, 
the photon spectrum is extremely suppressed. 

The qualitative behavior of the photon spectrum can be understood in the following way. 
If the photon energy is high enough, then double-photon pair creation is so efficient that 
this process dominates the cascade. However, once the photon energy becomes much smaller 
than 0{ml/T), this process is kinematically blocked. Numerically, this threshold is about 
mg/22T, as we mentioned. Then, photon- photon scattering dominates. However, since the 
scattering rate due to this process is proportional to E^, photon-photon scattering becomes 
unimportant in the limit — > 0. Therefore, for Ej <^ 0{ml/T), the remaining processes 
(pair creation in nuclei and inverse Compton scattering) are the most important. 

The crucial point is that the scattering rate for E^ ^ rn^/22T is much larger than that for 
E^ < ml/22T, since the number of targets in the former case is several orders of magnitude 
larger than in the latter. This is why the photon spectrum is extremely suppressed for 
-B-y ^ mg/22T. As a result, if the X particle decays in a thermal bath with temperature 
T > mg/22(5 (where Q is the binding energy of a nuclide) then photodissociation is not 
effective. 



B. Abundance of Light Elements with X 

Once the photon spectrum is formed, it induces the photodissociation of the light nuclei, 
which modifies the result of SBBN. This process is governed by the following Boltzmann 
equation: 

— ^N^ dE^(7N^-^N'{E^)f^{E^) 

SBBN N' 

^Y^riN" / dE^aN"^^N{E^)f'y{E^), (18) 

N" •' 

where is the number density of the nuclei A^, and [dnN/dt]sBBN denotes the SBBN 
contribution to the Boltzmann equation. To take account of the photodissociation processes, 
we modified the Kawano code [33], and calculated the abundances of the light elements. The 
photodissociation processes we included in our calculation are listed in Table II. 

The abundances of light nuclides will be functions of the lifetime of X (tx), the mass of 
X (mx), the abundance of X before electron-positron annihilation 

Yx = nx/n^, (19) 

and the baryon-to-photon ratio {rj). In our numerical BBN simulations, we found that 
the nuclide abundances depend only on the mass abundance mxYx, not on mx and Yx 
separately. In Figs. 5 - 9, we show the abundances of light nuclei in the mxYx vs. tx plane, 

at fixed 1]. 

We can understand the qualitative behaviors of the abundances in the following way. 
First of all, if the mass density of X is small enough, then the effects of X are negligible, 
and hence we reproduce the result of SBBN. Once the mass density gets larger, the SBBN 
results are modified. The effects of X strongly depend on tx, the lifetime of X. As we 



dn 



N 



dt 



+ SHun — 



dn 



N 



dt 
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mentioned in the previous section, photons with energy greater than ~ m1/22T participate 
in pair creation before they can induce photofission. Therefore, if the above threshold energy 
is smaller than the nuclear binding energy, then photodissociation is not effective. 

If Tx % 10"^ sec, then ml/22T < 2MeV at the decay time of X, and photodissociation is 
neghgible for all elements. In this case, the main effect of X is on the ^He abundances: if the 
abundance of X is large, its energy density speeds up the expansion rate of the universe, so 
the neutron freeze-out temperature becomes higher. As a result, ^He abundance is enhanced 
relative to SBBN. 

If 10^ sec <Tx ^ 10^ sec, then 2 MeV < ml/22T < 20 MeV. In this case, ^He remains 
intact, but D is effectively photodissociated through the process D + 7 — > p + n. When 
Tx ^ 10^ sec, ml/22T ^7.7 MeV (the binding energy of ^He), so ^He is dissociated for 
Tx ~ 10^ sec and large enough abundances mxYx ^ 10^^ GeV. However, D is even more 
fragile than ^He, so the ratio ^He/D actually increases relative to SBBN in this region, since 
it is dominated by D destruction. If the lifetime is long enough (rx ^ 10^ sec), ^He can also 
be destroyed effectively. In this case, the destruction of even a small fraction of the ^He can 
result in significant production of D and "'He, since the ^He abundance is originally much 
larger than that of D. This can be seen in Figs. 5 and 6: for Tx ^ 10^ sec and lO^^'' GeV 
^ mxYx ^ 10~^ GeV, the abundance of D changes drastically due to the photodissociation 
of '^He. Moreover, two-body decays of '^He into ^He or T (which decays into ^He) are preferred 
over the three-body decay ^He -I-7— >p-|-n-|-D, so the ^He/D ratio increases, relative to 
SBBN. If mxYx is large enough, all the light elements are destroyed efficiently, resulting in 
very small abundances. 

So far, we have discussed the theoretical calculation of the light-element abundances in 
a model with X decay. In the next section, we compare the theoretical calculations with 
observations, and derive constraints on the properties of X. 



C. Comparison with Observation 

Now, we compare the theoretical calculations with the observed abundances and show 
how we can constrain the model parameters. As we mentioned in Section II A, we have two 
^He values which are inferred from various observed data to be the primordial components. 
We will consider both cases and derive a constraint. (The statistical analysis we use to 
calculate the confidence level is explained in the appendix.) 

Let us start our discussion with the low ^He case {Y°^' = 0.234± (0.002)stat± (0.005),j^,t). 
Recalling that the low observed ^He value did not result in a good fit in the case of SBBN, 
we search the parameter space for regions of better fit than we can obtain with SBBN. 

In Fig. 10, we show the contours of the confidence level computed using four abundances 
(D/H, ^He, ^He/D, and ^Li/H), for some representative rj values (7710 = 2, 4, 5, 6), where 

7710 = 77 X 10^°. (20) 

The region of parameter space which is allowed at the 68% C.L. extends down to low 77 
(see Fig. 10a). Near 7710 = 2, deuterium is destroyed by an order of magnitude (without net 
destruction of "^He), so that the remaining deuterium agrees with the calculated low "^He. 
For 7/10 ~ 5, SBBN {rnxYx = 0) is allowed, and we have an approximate upper bound on 



10 



iTT-xYx (although for tx < 10^ sec, shghtly higher values of mxYx are allowed at low rj). 
For r]ic, ^ 6, no region is allowed at the 95% C.L., because 1] becomes to high to match even 
the observed D. We also plotted the regions excluded by the observational upper bounds on 
^Li/'^'Li. The shaded regions are i/e/?/? ^ 0.5, and the darker shaded regions are ye/vr ^ 1-3. 
Even if we adopt the stronger bound ye/y? ~ 0.5, our theoretical results are consistent with 
the observed ^Li value. 

In Fig. 11, we show the contours of the confidence levels for various lifetimes, tx = 
10^, 10^, 10^ sec. As the lifetime decreases, the background temperature at the time of decay 
increases, so the threshold energy of double-photon pair creation decreases. Then for a fixed 
mxXx", the number of photons contributing to D destruction decreases. Thus, for shorter 
lifetimes, we need larger mxYx in order to destroy sufficient amounts of D. The observed 
abundances prefer non- vanishing mxYx- 

In Fig. 12, we show contours of which have been projected along the rj axis into the 
Tx - nixYx plane. By projection, we mean taking the lowest C.L. value along the rj axis for 
a fixed point (tx, mxYx). The region above the sohd hke is excluded at the 95% C.L., while 
only the region within the dotted line is allowed at the 68% C.L. The 95% C.L. constraint 
for Tx ~ 10^ sec comes primarily from destruction of too much D; for tx ^ 10^ sec, it comes 
primarily from overproduction of ^He in ^He photofission. 

The lower mxYx region, i.e. mxYx ~ 10~^^ GeV, corresponds to SBBN, since there are 
not enough photons to affect the light-element abundances. It is notable that these regions 
are outside of the 68% C.L. This fact may suggest the existence of a long-lived massive 
particle X and may be regarded as a hint of physics beyond the standard model or standard 
big bang cosmology. 

For example, in Fig. 13 we show the predicted abundances of ^He, D, ^Li and ^Li adopting 
the model parameters tx — 10^ sec and mxYx = 5 x lO"^*^ GeV. The predicted abundances 
of ^He and ''Li are nearly the same as in SBBN. Only D is destroyed; its abundance decreases 
by about 80%. At low 7] ~ (1.7 — 2.3) x 10^^'^ in this model, the predicted abundances of 
these three elements agree with the observed values. It is interesting that the produced 
^Li abundance can be two orders of magnitude larger than the SBBN prediction in this 
parameter region. The origin of the observed ^Li abundance, ^Li/H ~ O(10~^^) is usually 
explained by cosmic ray spallation; however, our model demonstrates the possibility that 
^Li may have been produced by the photodissociation of ^Li at an early epoch. Our ^Li 
prediction is consistent with the upper bound Eq. (12). 

Although nixYx ^ 10^^° GeV is favored, it is worth noting that SBBN lies within the 
95% C.L. agreement between theory and observation. In Fig. 12, the 95% bound for tx ^ 10^ 
sec comes from the constraint that not much more than 90% of the deuterium should be 
destroyed; for tx ^ 10^ sec the constraint is that deuterium should not be produced from 
^He photofission. In Table III, we show the representative values of mxYx which correspond 
to 68% and 95% confidence levels respectively, for tx = 10'' — 10^ sec. 

Next, we would hke to discuss the high ^He case {Y°^' = 0.244 ± (0.002)5t„t± (0.005),^,*). 
Since the D abundance (6) and high value of ^He (10) both prefer a relatively high value 
of r], the SBBN prediction can be consistent with observation in this case. Therefore, we 
expect to be able to constrain the model parameters. 

For four representative r] values {rjiQ = 2, 4, 5, 6), we plot the contours of the confidence 
level in Fig. 14. In Fig. 2, we see that the SBBN calculations agree with the observed 



11 



abundances for mid-range values of the baryon-to-photon ratio (r/ ~ 5 x 10"''^'^). Thus, the 
approximate upper bound for rrixYx is plotted in Fig. 14c. (Again, for Tx < 10^ sec, slightly 
higher values of mxiOs: are allowed at low r}.) In Fig. 15, we show the C.L. plots for three 
typical lifetimes, tx — 10^, 10^, 10^ sec. This plot shows that SBBN works at better than 
the 68% C.L. for a range of hfetimes, but the non-standard scenarios with large mxYx and 
small 77 do not work as well as they did in the low Y°^^ case. Finally, we show the C.L. 
contours projected along the rj axis into the tx - mxYx plane (Fig. 16). Table IV gives the 
upper bounds on mxYx (GeV) which correspond to 68% and 95% C.L., for some typical 
values of the hfetime. 

Before we discuss additional constraints, let us comment on the case of high values of 
D/H, suggested by old observations [20-23]. Even though the high values seem less reliable, 
wc believe their possibility has not been completely ruled out. Therefore, it may be useful 
to comment on this case. The high value of D abundance [7/2 ^ O(10~^)] prefers a low value 
of 77, and hence it is completely consistent with the low value of Y°^^ in SBBN. Furthermore, 
if we adopt the relatively large error bar for y2 suggested by the observation, SBBN may 
also be consistent with high Y"^"^. Then, in this case also, we can obtain upper bounds on 
the mass density mxYx as a function of its lifetime. The upper bound behaves like the case 
of low D and low ^He shown in Fig. 16, and the upper bounds are within a factor of ten for 
most of lifetimes. 



We now mention additional constraints on our model. Since the cosmic microwave 
background radiation (CMBR) has been observed by COBE [48] to very closely follow 
a blackbody spectrum, one may be concerned about the constraint this gives on parti- 
cles with lifetime longer than ~ 10^ sec [49], which is when the double Compton process 
(7-1-6" # ^ + J + e~) freezes out [50].'' After this time, photon number is conserved, 
so photon injection from a radiatively decaying particle would cause the spectrum of the 
CMBR to become a Bose-Einstein distribution with a finite chemical potential fi. COBE [48] 
observations give us the constraint \fi\ ^ 9.0 x 10~^. For small fi, the ratio of the injected 
to total photon energy density is given by 6p^/p^ ~ 0.71/i.. Thus, we have the constraint 



Note that the CMBR constraint is not as strong as the bounds we have obtained from 
BBN. In particular, ^He/D gives us our strongest constraint for lifetimes longer than 10^ 
sec, because ^He photofission overproduces ^He [26]. 

In this paper, we have considered only radiative decays; i.e., decays to photons and in- 
visible particles. If X decayed to charged leptons, the effects would be similar to those of the 



^This constraint applies only to particles with lifetime shorter than ~ 4 x 10^° sec, which corre- 
sponds to the decoupling time of Compton/inverse Compton scattering. After this time, injected 
photons do not thermalize with the CMBR. 



D. Additional Constraints 




(21) 
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decay to photons, because charged leptons also generate electromagnetic cascades, resulting 
in many soft photons. On the other hand, if X decayed only to neutrinos, the constraints 
would become much weaker. In the minimal supersymmetric standard model (MSSM), the 
X particle would decay to neutrinos and sneutrinos. The emitted neutrinos would scatter off 
of background neutrinos, producing electron-positron pairs, which would trigger an electro- 
magnetic cascade. Because the interaction between the emitted and background neutrinos is 
weak, the destruction of the light elements does not occur very efficiently [51]. In contrast, if 
X decayed to hadrons, we expect that our bounds would tighten, because hadronic showers 
could be a significant source of D, ^He, ^Li, ^Li, and ^Be [6]. In fact, even though we have 
assumed that X decays only to photons, these photons may convert to hadrons. Thus, the 
branching ratio to hadrons is at least of order 1 %, if kinematically allowed [7]. Since we 
have neglected this effect, our photodissociation bounds are conservative. 

IV. MODELS 

So far, we have discussed general constraints from BBN on radiatively decaying particles. 
In the minimal standard model, there is no such particle. However, some extensions of the 
standard model naturally result in such exotic particles, and the light-element abundances 
may be significantly affected in these cases. In this section, we present several examples of 
such radiatively decaying particles, and discuss the constraints. 

Our first example is the gravitino ip, which appears in all the supergravity models. The 
gravitino is the superpartner of the graviton, and its interactions are suppressed by inverse 
powers of the reduced Planck scale ~ 2.4 x 10^^ GeV. Because of this suppression, the 
lifetime of the gravitino is very long. Assuming that the gravitino's dominant decay mode 
is to a photon and its superpartner (the photino), the gravitino's lifetime is given by 

T3/2 ~ 4 X 10^ sec X (m3/2/l TeV)-^ (22) 

where 777.3/2 is the gravitino mass. Notice that the gravitino mass is 0(100 GeV — 1 TeV) 
in a model with gravity-mediated supersymmetry (SUSY) breaking, resulting in a lifetime 
which may affect BBN. 

If the gravitino is thermally produced in the early universe, and decays without being 
diluted, it completely spoils the success of SBBN. Usually, we solve this problem by introduc- 
ing inflation, which dilutes away the primordial gravitinos. However, even with inflation, 
gravitinos are produced through scattering processes of thermal particles after reheating. 
The abundance I3/2 = 77,3/2/n^ of the gravitino depends on the reheating temperature Tr, 
and is given by [4] 

^3/2 ^ 3 X 10-" X (TR/10^°GeV). (23) 

Therefore, if the reheating temperature is too high, then gravitinos are overproduced, and 
too many light nuclei are photodissociated. 

We can transform our constraints on (rx, 777x^x) into constraints on (7773/2, Tr). In 
particular, we use the projected 95% C.L. boundaries from Figs. 12 and 16. For several 
values of the gravitino mass, we read off the most conservative upper bound on the reheating 
temperature from Fig. 17, and the results are given by 
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m3/2 = 100 GeV (r3/2 ~ 4 x 10^ sec) : Tr < 2 x 10^ GeV, 
m3/2 = 1 TeV (t3/2 ~ 4 x 10^ sec) : Tr < 6 x 10^ GeV, 
m3/2 = 3 TeV (T3/2 ~ 1 x 10^ sec) : Tr < 2 x 10^^ GeV. 

If the gravitino is heavy enough (m3/2 ^ 5 TeV), then its hfetime is too short to destroy 
even D. In this case, our only constraint is from the overproduction of ^He. If the gravitino 
mass is hghtcr, then the hfetime is long enough to destroy D, ^He, or even ^He. In this case, 
our constraint on the reheating temperature is more severe. 

Another example of our decaying particle is the lightest superparticle in the MSSM 
sector, if it is heavier than the gravitino. In particular, if the lightest neutralino is the 
hghtest superparticle in the MSSM sector, then it can be a source of high-energy photons, 
since it wih decay into a photon and a gravitino. In this case, we may use BBN to constrain 
the MSSM. 

The abundance of the lightest neutralino is determined when it freezes out of the thermal 
bath. The abundance is a function of the masses of the superparticles, and it becomes larger 
as the superparticles get heavier. Thus, the upper bound on nixYx can be translated into 
an upper bound on the mass scale of the superparticles. 

In order to investigate this scenario, we consider the simplest case where the lightest 
neutralino is (almost) purely bino B. In this case, the lightest neutralino pair-annihilates 
through squark and slepton exchange. In particular, if the right-handed sleptons are the 
lightest sfermions, then the dominant annihilation is B + B — > l'^ + l" . The annihilation 
cross section though this process is given by [52] 

9, 9x f "^1 2m| 2m% 1 ^ ^ 

where (v^) is the average velocity squared of bino, and we added the contributions from 
all three generations by assuming the right-handed sleptons arc degenerate.^ With this 
annihilation cross section, the Boltzmann equation for the number density of binos is given 
by 

hs + SHn^ = -2((7T;,e,)(4 " (riff), (25) 

where is the equilibrium number density of binos. The factor 2 is present because two 
binos annihilate into leptons in one collision. We solved this equation and obtained the mass 
density of the bino as a function of the bino mass and the right-handed slepton mass. (For 
details, see e.g. Ref. [53]). Numerically, for = 100 GeV, mx^x ranges from ~ 10~^ GeV 
to ~ 10~^ GeV as we vary mj from 100 GeV to 1 TeV. If mxYx is in this range, BBN 
is significantly affected unless the lifetime of the bino is shorter than 10^ - 10^ sec (see 
Tables 111 - IV). The lifetime of the bino is given by 



^If the bino is heavier than the top quark, then the s-wave contribution annihilating into top 
quarks becomes important. In this paper, we do not consider this case. 
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Notice that the Ufetime becomes shorter as the gravitino mass decreases; hence, too much 
D and ''Li are destroyed if the gravitino mass is too large. Therefore, we can convert the 
constraints given in Figs. 12 and 16 into upper bounds on the gravitino mass. Since the 
abundance of the bino is an increasing function of the slepton mass m^-^, the upper bound on 
the gravitino mass is more severe for larger slepton masses. For example, for = 100 GeV, 
the upper bound on the gravitino mass is shown in Fig. 18. At some representative values 
of the slepton mass, the constraint is given by 

mj^ = 100 GeV : m3/2 ^ 1 GeV, 
rrii^ = 300 GeV : mg/s < 700 MeV, 
mj^ = 1 TeV : rng/a < 400 MeV. 

As expected, for a larger value of the slepton mass, the primordial abundance of the bino 
gets larger, and the upper bound on the gravitino mass becomes smaller. 

Another interesting source of high-energy photons is a modulus field 0. Such fields are 
predicted in string-inspired supergravity theories. A modulus field acquires mass from SUSY 
breaking, so we estimate its mass m<^ to be of the same order as the gravitino mass (see for 
example [54]). 

In the early universe, the mass of the modulus field is negligible compared to the expan- 
sion rate of the universe, so the modulus field may sit far from the minimum of its potential. 
Since the only scale parameter in supergravity is the Planck scale M*, the initial amplitude 
4>o is naively expected to be of 0(M*). However, this initial amplitude is too large; it leads 
to well-known problems such as matter domination and distortion of the CMBR. Here, we 
regard 0o as a free parameter and derive an upper bound on it. 

Once the expansion rate becomes smaller than the mass of the modulus field, the modulus 
field starts oscillating. During this period, the energy density of is proportional to 
(where R is the scale factor); hence, its energy density behaves like that of non-relativistic 
matter. The modulus eventually decays, when the expansion rate becomes comparable to 
its decay rate. Without entropy production from another source, the modulus density at 
the decay time is approximately 

m^Y^ = ^ - 5 X 10^° GeV x (m^/1 TeVy/\(Po/M,f , (27) 

where is the energy density of the modulus field. As in our other models, we can convert 
our constraints on {Tx,mxYx) (Figs. 12, and 16) into constraints on (m<^,0o)- Using the 
most conservative of these constraints, we still obtain very stringent bounds on the initial 
amplitude of the modulus field 

= 100 GeV (r^ ~ 4 x 10^ sec) 
= 1 TeV (r^ ~ 4 X 10^ sec) 
= 3 TeV (r^ ~ 1 x 10^ sec) 



487r 



m| cos^ 6'w 



^0 ~ -I- X iu «jev , 
00 ^ 5 X 10^ GeV, 
6n<9x 10^ GeV. 
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Clearly, our upper bound from BBN rules out our naive expectation that 0o ~ M*. It is 
important to notice that (conventional) inflation cannot solve this difliculty by diluting the 
coherent mode of the modulus field. This is because the expansion rate of the universe 
is usually much larger than the mass of the modulus field, and hence the modulus field 
does not start oscillation. One attractive solution is a thermal infiation model proposed 
by Lyth and Stewart [55]. In the thermal inflation model, a mini-inflation of about ~ 10 
e-folds reduces the modulus density. Even if thermal inflation occurs, there may remain a 
significant modulus energy density, which decays to high-energy photons. Thus, BBN gives 
a stringent constraint on the thermal infiation model. 

V. DISCUSSION AND CONCLUSIONS 

We have discussed the photodissociation of light elements due to the radiative decay of 
a massive particle, and we have shown how we can constrain our model parameters from 
the observed light-element abundances. We adopted both low and high '^He values in this 
paper, and we obtained constraints on the properties of the radiatively decaying particle in 
each case. 

When we adopt the low ^He value, we find that a non-vanishing amount of such a long- 
lived, massive particle is preferred: rrixYx ^ 10^^'^GeV and lO^sec ^ Tx ^ lO^sec. On 
the other hand, consistency with the observations imposes upper bounds on mxXx in each 
cases. 

We have also studied the photodissociation of ^Li and ^Li in this paper. These processes 
do not affect the D, ^He, and ^He abundances, because ^Li and ^Li are many orders of 
magnitude less abundant than D,"^He, and '^He. When we examine the region of parameter 
space where the predicted abundances agree well with the observed ^Li and the low "^He 
observations, we find that the produced ®Li/H may be of order 10~^^, which is two orders 
of magnitude larger than the prediction of SBBN (see Figs. 7 and 13). The predicted ^Li is 
consistent with the observed upper bound Eq. (12) throughout the region of parameter space 
we are interested in. Although presently it is believed that the observed ®Li abundance is 
produced by spallation, our model suggests another origin: the observed ®Li may be produced 
by the photodissociation of ^Li. 

We have also discussed candidates for our radiatively decaying particle. Our first exam- 
ple is the gravitino. In this case, we can constrain the reheating temperature after infiation, 
because it determines the abundance of the gravitino. We obtained the stringent bounds 
Tr < lO^GeV — lO^GeV for 100 GeV ^ "023/2 ^ 1 TeV. Our second example is the lightest 
neutralino which is heavier than the gravitino. When the neutralino is the lightest super- 
particle in the MSSM sector, it can decay into a photon and a gravitino. If we assume the 
lightest neutralino is pure bino, and its mass is about 100 GeV, then the relic number density 
of binos is related to the right-handed slepton mass, because they annihilate mainly through 
right-handed slepton exchange. For this case, we obtained an upper bound on the gravitino 
mass: 777.3/2 ^ 400 MeV — 1 GeV for 100 GeV ^ 777^-^ ^ 1 TeV. Our third example is a mod- 
ulus field. We obtained a severe constraint on its initial amphtude, 0o ~ 10^ GeV — 10^ GeV 
for 100 GeV ^ 7773/2 ^ 1 TeV. This bound is well below the Planck scale, so it suggests the 
need for a dilution mechanism, such as thermal infiation. 



16 



ACKNOWLEDGEMENT 



This work was supported by the Director, Office of Energy Research, Office of Basic 
Energy Services, of the U.S. Department of Energy under Contract DE-AC03-76SF00098. 
K.K. is supported by JSPS Research Fellowship for Young Scientists. 



APPENDIX 

In this appendix, we explain how we answer the question, "How well does our simulation 
of BBN agree with the observed light-element abundances?" To be more precise, we rephrase 
the question as, "At what confidence level is our simulation of BBN excluded by the observed 
light-element abundances?" 

From our Monte-Carlo BBN simulation, we obtain the theoretical probability density 
function (p.d.f.) p^^j{yl^,Y^^,log'^(^yf^) of our simulated light-element abundances yl^,Y^^, 
and log^Q?/*'^. Wc find that P247(?/2^> ^^'^j ^^Sio well- approximated by a multivariate 

Gaussian distribution function: 

pIM, y'\ logio yf) = "''(yf) X pf "-(y*'^) x p^^-'^iog^^ yf), (Ai) 

where 

1 — ^ 



Gauss/ ^ ^\ ^ J_ exp 



2\ a 



(A2) 



Note that p^^^iy^, Y^^, log^Q yf^) depends upon the parameters p = (77, ...) of our theory. 
(The ellipses refer to parameters in non-standard BBN, e.g., ruxYx, Tx.) In particular, the 
means and standard deviations of p^^j{yl^, Y^^, logi^y^^) arc functions of p. 

For BBN with a radiatively decaying particle, we also consider the ratio r*'' = yl^/y^^. 
Our Monte-Carlo BBN simulation allows us to find the p.d.f. p||(r*'*). We approximate 
by a Gaussian, and neglect the correlation between r and y^, Y^^, logio Ut^- This assumption 
(which is justified in work to be published by one of the authors [24]), allows us to write 

p%M\ logio P) = P'^{r'\ p) X p*'47«, logio ?/f ; P)- (A3) 

We want to compare these theoretical calculations to the observed light-element abun- 
dances y^'^^Y"^^, and Xogi^yf'^. Since the observations of the light-element abundances are 
independent, we can factor the p.d.f. p'^l^{yf'^,Y"^%\og^Qyf'^) as 

pt,{yt. y'', logio yr'n = Ptiyt) X Pt{y"'') X (logio yf')- (A4) 

We assume Gaussian p.d.f. 's for 7/2''^, Y°'^^, and log^g y?^^- We use the mean abundances and 
standard deviations given in Equations (6), (9), (10), and (11). Since we have two discordant 
values of ^Hc, we considered two cases, i.e., high and low values of ''He abundances. 

^He is more complicated. Aside from the trivial positivity requirement, we only have an 
upper bound on r°^* (the primordial ^He/D, as deduced from solar-system observations); 
namely, r°^^ < Tq. Because of this, we choose for the p.d.f. 
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^obs („obs 



r , r°^* < 



N , < r"''' < 



^ j.obs 



(A5) 



where the normaUzation factor is = l/{r^ + 0-^0^71/2), and r^, 0-^0 are given in Eq. (7). 

Since the observations of the hght-element abundances are independent, we can write 
the total observational p.d.f. for BBN+X as 

P2U7(y2'', y'', r"'', logio Vv"') = P2f{r'"') x P^%{yt, y'', logio Vv"')- (A6) 
To simplify the notation, we write 

= (^yl\Y'\r'Mog,,y',') (A7) 

b*" = (l/*^y*^loglol/f) (A9) 

h"'^ = {yt,Y'"^log,,y",'^). (AlO) 
Then the quantities Aa = a*^ — a°^^ have a p.d.f. given by: 

p^47(Ab) = J dh"^' P247(b°'") J dh*'' p*^7(b*'^)5(Ab - (b*'* - h"^')) 

^ jdh pfMpfi,{h - Ab), [SBBN] (All) 

= I da p*2347(a)p^347(a - Aa), [BBN+X] (A12) 

where we have suppressed the dependence of ^^^(Ab), P247(b*^)) P2347(Aa), and P2347(a*'*) 
on the theory parameters p. The integral in Eq. (All) is simply a Gaussian: 



p^47(Ab) = n /?r- 



i=2,4,7 V27r(7i 



1 / Abi - A6,- 



2 V a\ 



.2 



(A13) 



where Abi = bf' — bf^, af = {aj'^)^ + {af'^)^, and i runs over y2, Y and logioy7. To evaluate 
Eq. (A12), we note that 

P2347(a°'1 = P23{r"'') X plfri^^'n (AM) 

P^'347(a*'^) = P23(r*'^) X p*'47(b*'^), (A15) 

where P247 P247 Gaussian. The integral then factors as 

p^347(Aa) = I ^^rp*^3(0P23'(^ -Ar)x j dh pt,{h)pfi,{h - Ab) (A16) 
= p^3(Ar)xp^4,(Ab). (A17) 

The first integral can be evaluated as 
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(A18) 



-oo 



dr , — — exp 



+ 

N 
~2 



dr 



N 



1 /r + Ar — r*'* 

2 l~ a, 



th 



exp — - 



V + Ar — r 



erf I 

2 



Vq + Ar — r 
V2^ 



— erf 



+ 



^0 



rO 



(A19) 



th 



1 — erf 



^0 



exp 



1 /Ar-r^^ + f^' 

2 I 



(A20) 



where a,"' = (af + (a,©)-^, f/cf,^ = /^/(a.o)^ - (Ar - r^f^) / {a',^)^ , and = (af )2 + 

Our question can now be rephrased as, "At what confidence level is Aa = excluded?" 
To answer this, we need to consider the region 5" in abundance space where the value of the 
p.d.f. is higher than 



P 



r p^47(Ab = 0;p) [SBBN] 
\p2'^347(Aa=0;p) [BBN+X] 



Mathematically phrased, 

C.L.(p) = / rf(Ab)p^47(Ab;p), [SBBN] 



where 



C.L.(p) = / d(Aa)p|47(Aa;p), [BBN+X] 



S = {Ah:p^,,{Ah;p)>p}, [SBBN], 
5 = {Aa : p^347(Aa; p) > p}, [BBN+X]. 



(A21) 

(A22) 
(A23) 



(A24) 
(A25) 



We use the CL. to constrain various scenarios of BBN. 

In the SBBN case, the integral is Gaussian and is easily expressed in terms of — 
-21og[(27r)3/V2(T4(T7p] (See Eqs. (15) and (17)): 



C.L.(p) 



TT 



exp 



+ erf 



[SBBN]. 



(A26) 



To evaluate the BBN+X integral Eq. (A23), we separate the Gaussian variables Ab from 
the non-Gaussian variable Ar, using the decomposition in Eq. (A16): 



C.L.(p)=/ d{Ar)pUAr) [ d{Ah)pU^h), 

JS23 JS2i7{Ar) 



(A27) 



where 5'23 = the set of Ar such that p^(Ar) > p/max(p247); and 5'247(Ar) = the set of Ab 
such that p^7(Ab) > p/p^(Ar). We can easily evaluate the Gaussian integral. Again, the 
result is conveniently expressed in terms of x^. 
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( 



x(Ar) 



) 



(A28) 



where x(A^) = y — 2 log[(27r)3/2o'20"40"7p/j923(Ar)]. We then evaluate Eq. (A28) numerically. 



Our confidence level is calculated for four degrees of freedom Aaj (three, in the case 
of SBBN). It denotes our certainty that a given point p in the parameter space of the 
theory is excluded by the observed abundances. In order to compare our theory with a 
late-decaying particle (three parameters p: rxjirixYx, and rj) to a theory with a different 
number of parameters {e.g., only one in SBBN), one would want to use a variable in these 
parameters. This transformation would be possible if the abundances were linear in the 
theory parameters p. In this case, we could integrate out a theory parameter such as rj and 
set a C.L. exclusion limit (with a reduced number of degrees of freedom) on the remaining 
parameters. However, the turn out to be highly non-linear functions of p, so integrating 
out a theory parameter turns out to have little meaning. Instead, we project out various 
theory parameters (as explained in Section III C) to present our results as graphs. 
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TABLES 




iV^ (95 


% C.L.) 


7] X IQio (95 % C.L.) 


Low "^He 


2.1; 


fl.O 
-0.8 


^•'-0.8 


High ^He 


2.8: 


fl.O 
-1.0 


5 

"^•u-0.8 



TABLE L Observational constraints on 77 and ATj, in SBBN 





Photofission Reactions 


Icr Uncertainty 


Threshold Energy 


Ref. 


1. 


D + 7 — > p + n 


6% 


2.2 MeV 


[40] 


2. 


T + 7^ n + D 


14% 


6.3 MeV 


[41,42] 


3. 


T + 7 ^ p + 2n 


7% 


8.5 MeV 


[42] 


4. 


^He + 7 ^ p + D 


10% 


5.5 MeV 


[43] 


5. 


^He + 7 — n + 2p 


15% 


7.7 MeV 


[43] 


6. 


^He + 7 ^ p + T 


4% 


19.8 MeV 


[43] 


7. 


^He + 7 ^ n + -"^Hc 


5% 


20.6 MeV 


[44,45] 


8. 


^He + 7^p + n + D 


14% 


26.1 MeV 


[46] 


9. 


^Li + 7 ^ anything 


4% 


5.7 MeV 


[47] 


10. 


^Li + 7 — ^ 2n + anything 


9% 


10.9 MeV 


[47] 


11. 


'''Li + 7 ^ n + ^Li 


4% 


7.2 MeV 


[47] 


12. 


''Li + 7 ^ ^He + anything 


9% 


2.5 MeV 


[47] 


13. 


''Be + 7 ^ p + ""Li 








14. 


^Be + 7 ^ anything except ^Li 









TABLE II. Photodissociation processes, and the 1-a uncertainty in the cross sections. Since 
there is no experimental data on photodissociation of ^Be, we assume in this paper that the rate 
for Reaction 13 is the same as for Reaction 11, and the rate for Reaction 14 is the sum of the rates 
for Reactions 10 and 12. 
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Tx = 10^ sec 


10^ sec 


106 sec 


10^ sec 


10^ sec 10^ sec 


95% C.L. 6x10-6 


9x10-9 


1x10-*^ 


3x10-12 


4x10-13 3x10-13 


68% C.L. (6-60)xl0-'^ 


(8 - 80) X 10-1° 


(1 - 9) X 10-1° 







TABLE III. Upper or (upper - lower) bound on mx^^ in units of GeV for the case of low 
^He. Note that the C.L. is for four degrees of freedom, and r) is varied to give the extreme values 
for mxYx- 





Tx = 


10^ sec 


10^ sec 


106 sec 


10^ sec 


10^ sec 


109 sec 


95% C.L. 


7x 


10-6 


7x 10-9 


8x 10-1° 


6x 10-12 


7x 10-13 


5x 10-13 


68% C.L. 


5x 


10-6 


5x 10-9 


6x 10-1° 


3x 10-12 


5x 10-13 


3x 10-13 



TABLE IV. Same as Table III, except for high ^He. 
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FIG. 1. SBBN prediction of the abundances of the hght elements. The sohd hnes are the central 
values of the predictions, and the dotted lines represents the 1-a uncertainties. The boxes denote 
the 1-a observational constraints. 
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FIG. 2. &s a function of rj, for SBBN with three degrees of freedom (r/, tx, w-x^)- We show 
our results for both of the ^He abundances deduced from observation: low '^He (dashed), high ^He 
(solid). 
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FIG. 3. Figure 3: C.L. for BBN as a function of r/ and N^, with (a) low value of Y , and (b) 
high value of Y . The filled square denotes the best-fit point. 



27 




10-3 10-2 10-1 1 101 102 103 



Energy (GeV) 
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FIG. 5. The abundance of D/H in the mxYx vs. tx plane with (a) = 2 x 10 (b) 
r] = 4x 10-1°, (c) r/ = 5 x 10"^°, and (d) r/ = 6 x 10"^°. 
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FIG. 6. The mass fraction of ^He, for the same theory parameters as in Fig. 5. 
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FIG. 7. The abundance of ^Li/H, for the same theory parameters as in Fig. 5. 
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FIG. 8. The abundance of ^He/D, for the same theory parameters as in Fig. 5. 
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FIG. 9. The abundance of '^Li/H, for the same theory parameters as in Fig. 5. 
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FIG. 10. C.L. in the mx^x vs. tx plane, for low value of Y . We take (a) ?7 = 2 x 10^^*^, (b) 
r/ = 4 X 10^^°, (c) r/ = 5 X 10-i^°, and (d) 77 = 6 x IQ-i^^. The shaded regions are y^jy-j > 0.5, and 
the darker shaded regions are ye/y? ^ 1-3. 
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FIG. 11. C.L. in the rj vs. mxYx plane for various values of tx, for low value of Y. 
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FIG. 12. Contours of C.L. projected on r] axis for low value of Y. 
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FIG. 13. Predicted light-element abundances ^He, D, ^Li and ^Li at tx = 10^ sec and 

mxYx = 5 X 10^"^*^ GcV. The contours which are favored by observation are plotted, adopting 
the low '^Hc and low D values. The dotted line denotes the 95% C.L. and the shaded region de- 
notes the 68% C.L.. The predicted ^Li abundance is two orders of magnitude larger than the case 
of SBBN. 
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FIG. 14. Same as Fig. 10, except for high value of Y. 
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FIG. 15. Same as Fig. 11, except for high vakie of Y. 
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FIG. 16. Same as Fig. 12, except for high value of Y. The region above the solid line is excluded 
at the 95% C.L., while the region above the dotted line is excluded at the 95% C.L. 
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FIG. 17. Contours of 95% C.L., yielding an upper bound on the reheating temperature, as a 
function of the gravitino mass. These are the results of the two-photon emission by the decay of a 
gravitino. If we consider the one-photon emission such as V'/^ — ^ 7 + 7; the upper bounds become 
milder by a factor of two. 
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FIG. 18. Contours of 95% C.L., yielding an upper bound on the gravitino mass, as a function 
of the right-handed slepton mass. 
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